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Abstract 

Detecting "B-mode" (i.e., divergence free) polarization in the Cosmic Microwave Background (CMB) 
would open a new window on the very early Universe. However, the polarized microwave sky is 
dominated by polarized Galactic dust and synchrotron emissions, which may hinder our ability to test 
inflationary predictions. In this paper, we report on our knowledge of these "Galactic foregrounds," as 
well as on how a CMB satellite mission aiming at detecting a primordial B-mode signal ( "CMBPol" ) will 
contribute to improving it. We review the observational and analysis techniques used to constrain the 
structure of the Galactic magnetic field, whose presence is responsible for the polarization of Galactic 
emissions. Although our current understanding of the magnetized interstellar medium is somewhat 
limited, dramatic improvements in our knowledge of its properties are expected by the time CMBPol 
flies. Thanks to high resolution and high sensitivity instruments observing the whole sky at frequencies 
between 30 GHz and 850 GHz, CMBPol will not only improve this picture by observing the synchrotron 
emission from our galaxy, but also help constrain dust models. Polarized emission from interstellar dust 
indeed dominates over any other signal in CMBPol's highest frequency channels. Observations at these 
wavelengths, combined with ground-based studies of starlight polarization, will therefore enable us to 
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improve our understanding of dust properties and of the mechanism(s) responsible for the alignment 
of dust grains with the Galactic magnetic field. CMBPol will also shed new light on observations that 
are presently not well understood. Morphological studies of anomalous dust and synchrotron emissions 
will indeed constrain their natures and properties, while searching for fluctuations in the emission from 
heliospheric dust will test our understanding of the circumheliospheric interstellar medium. Finally, 
acquiring more information on the properties of extra-Galactic sources will be necessary in order to 
maximize the cosmological constraints extracted from CMBPol's observations of CMB lensing. 
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1 Introduction 



Over the course of the last two decades, studies of the Cosmic Microwave Background (CMB) 
have revolutionized our understanding of the composition, structure and evolution of the Uni- 
verse. After sixteen years of observations of the CMB temperature anisotropies first detected 



by the COBE satellite in the early 1990's (Smoot et al. 1992), we now have a standard cosmo- 



logical model fully described by a handful of parameters, and in spectacular agreement with 



a broad range of independent astrophysical observations (Komatsu et al. , 20081. 



The cornerstone of this achievement is doubtless inflation. This paradigm, initially de- 



signed to explain the observed flatness and isotropy of the Universe (Guth 19811, provides a 



natural explanation for several of its other properties, including the origin of the primordial 
density fluctuations that led to structure formation. A generic prediction of the most natural 
inflationary models is the existence of primordial gravitational waves that should have left 
their imprint on the CMB anisotropies in the form of a divergence-free polarization pattern 
known as "B-mode polarization." However, this signal is expected to have a very low ampli- 
tude, typically more than three to four orders of magnitude smaller than the measured CMB 
temperature fluctuations (for a review of inflationary predictions, see the companion paper by 



Baumann et al. , 2008a I . Detecting primordial B-modes is therefore a challenging enterprise 
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made only more difficult by polarized emissions from our galaxy and others. Among these so- 
called "polarized foregrounds" are polarized synchrotron and dust emissions from the Galaxy, 
which dominate the expected primordial B-mode signal over the whole sky at all frequencies 
of interest (see, e.g., the discussion in Page et al. 2007 1 . Thus, being able to properly handle 
them is central to being able to detect low amplitude B-modes in the CMB. 

In this document, we do not address the question of how to look behind this Galactic (and 
extra-Galactic) screen, but focus on our knowledge of the foregrounds themselves. Our ability 
to use the data provided by a satellite mission sensitive enough to detect B-modes of inflation- 
ary origin (hereafter, CMBPol) is indeed limited by our understanding of foregrounds (see the 
companion paper by Dunkley et al. 2008a). We therefore review what we know now, how we 
know it, and how we hope to improve it by the time CMBPol flies. Moreover, since Galactic 
foregrounds are the dominant signal CMBPol will observe, we can expect our understanding 
of these emissions and their underlying physical processes to significantly improve through 
an analysis of the CMBPol data. Although our current knowledge of foregrounds makes it 
difficult to quantify many of these possible improvements, we review what we think are the 
main ways in which CMBPol will help us explore the magnetized Galactic interstellar medium. 



Four companion papers complement the contents of this document: Baumann et al. (2008a 



review the inflationary science case for CMBPol; Smith et al. (20081 report on CMB lensing; 



Zaldarriaga et al. (20081 discuss reionization; and Dunkley et al. (2008a) provide a detailed 



analysis of how and how well a primordial B-mode signal can be separated from the polarized 



emissions from the Galaxy. An executive summary of these documents is also available ( Bau- 



mann et al. 2008b). For the purposes of this paper and all companion documents, we assume 



that CMBPol will have specifications close to those provided in Tab. [T] which fall between the 
specifications of the two EPIC designs proposed in the Bock et al.| ( |2008[ ) report. 



Table 1: Possible CMBPol'' Specifications 



Frequency (GHz) 


30 


45 


70 


100 


150 


220 


340 


500 


850 


Resolution (FWHM in arcmin) 


26 


17 


11 


8 


5 


3.5 


2.3 


1.5 


0.9 


Sensitivity'^ (nK CMB) 


160 


69 


35 


27 


26 


40 


180 


150^= 


400^^ 



^Expected mission lifetime of 4 years 
''CMB temperature fluctuations detectable at la in a 2° x 2° pixel 
■^Could be improved by a factor of 5 by cooling the telescope to 4 K 
''Could be improved by a factor of 40 by cooling the telescope to 4 K 



This document is organized as follows. In Sec. [2| we review the observational techniques 
used to probe the Galactic magnetic field (Sec. |2.1[ ), as well as what they teach us about its 
main properties (Sec. 2.2), while Sec. [3] presents an analysis technique aimed at improving 
our constraints on the magnetized interstellar medium. Sec. |4] discusses the progress one can 
expect from upcoming radio (Sec. 4.1) and CMB (Sec. [43] ) observations, including those by 
CMBPol. Our knowledge of the polarized emission from interstellar dust is then summarized 
in Sec. |5] and followed by a description of the physical mechanism responsible for it in Sec. |6j 
Heliospheric dust emission, and "anomalous emissions" detected in the WMAP data, are the 
topics of, respectively. Sec. [7] and Sec. [8j Finally, Sec. |9] gives an overview of the foreground 
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situation in the context of CMB lensing studies, and our conclusions are presented in Sec. 
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2 Probing the Galactic Magnetic Field 



The idea that the Galaxy has a magnetic field, and that it may play an important role 



in Galactic physics, dates back to Fermi (1949 1 . He proposed that cosmic rays might be 



generated outside of the solar system as a result of energetic particles colliding with moving 
irregularities in this field. In his view, the magnetic field would not only be a cosmic ray 
generator, but also a containment factor preventing the rays from escaping the Galaxy. 

It is now believed that magnetic fields and cosmic rays indeed contribute to the vertical 
support of the gas in the Galaxy, with the magnetic field energy density dominating the tur- 



bulent and thermal gas energy densities at high Galactic latitudes (Boulares Sz Cox, 19901 



However, astronomical magnetic fields also play important roles in many other processes af- 
fecting the structure of the Galaxy. In star formation, they can inhibit the gravitational 
collapse of interstellar clouds (primary star formation regions) and help removing prestellar 



angular momentum (Zweibel & Heiles, 1997), thereby affecting the distribution of stars in the 



Galaxy. It is also believed that magnetic fields influence galaxy formation and evolution by 



causing large density fluctuations which result in structures within galaxies (Kim et al. 19961. 



Yet despite the Galactic magnetic field's importance, many mysteries remain about its 
generation and evolution. A dynamo mechanism is believed to be necessary to increase galactic 
magnetic field strengths from that of a tiny primordial seed field to the observed values (see, 
e.g.. Beck et al. 19961, but neither the exact dynamo mechanism nor the generation and 



strength of the seed field are significantly constrained by available observations. External 



galaxies appear to have a wide range of magnetic field characteristics (Beck, 2001 ), indicating 



that different dynamos must be acting in different galaxies. Furthermore, in some cases, rapid 



field amplification seems to be required (Gaensler et al. 2005). Before we can hope to fully 



explore these aspects of the field, we must first understand the field's general topology. 

While observations of polarized emission from interstellar dust grains, polarization of 
starlight, and Zeeman splitting can provide in depth probing of the field in localized regions, 
most of the knowledge we have about the magnetic field across the bulk of the Galaxy comes 
from radio observations of linearly-polarized radiation. At frequencies above a few GHz, maps 
of radio polarization angles directly yield the direction of the magnetic field component in the 
plane of the sky, while the total emission indicates its strength. At lower frequencies, Faraday 
rotation of polarized emission traveling through the interstellar medium (ISM) can be used to 
probe the magnetic field component parallel to the line-of-sight. We discuss these techniques 



in Sec. 2.1 before detailing the insights gained from these observations in Sec. 2.2 



2.1 Observational Techniques 
2.1.1 Radio Synchrotron Radiation 

Relativistic electrons moving in the Galactic magnetic field emit synchrotron radiation. Opti- 
cally thin synchrotron emission, eg, at frequency v depends on the component of the uniform 
magnetic field perpendicular to the line-of-sight, B±, according to eg oc zy~(p^i)/2 



(Rybicki & Lightman 1979). The exponent p describes the power law energy spectrum of the 
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ensemble of relativistic electrons through the relation N(E) oc E~^. This emission therefore 
provides information about the slope of the electron energy spectrum, as well as an estimate 
of the strength of the magnetic field if assumptions are made about the volume of the source 
and the energy density of cosmic rays. The minimum combined energy density of cosmic ray 
particles and magnetic field necessary to reproduce the emission from an observed source is 
similar to the equipartition energy density between magnetic field and cosmic rays. Magnetic 
field estimates from the brightness of synchrotron emission usually assume this minimum 
energy condition, or equipartition, though often without clear observational evidence (through 
7- or X-ray data) that it applies (see, e.g., the discussion in Beck &; Krause 20051. 

Synchrotron emission from a region with a uniform magnetic field and an electron distri- 
bution with p = 3 (respectively, p = 2) has a theoretical polarization of 75% (respectively, 
60%), with the plane of polarization perpendicular to the direction of B± in the plane of the 
sky ( [Rybicki Lightman 1979). Polarization of synchrotron emission therefore also gives in- 
formation on the magnetic field component perpendicular to the line-of- sight. In practice, the 
observed emission is integrated over large regions in space with a complicated magnetic field 
structure and an ionized plasma present in addition to the relativistic electrons. Consequently, 
the integrated emission is usually much less polarized than the theoretical limit. 

Although maps of the synchrotron emission can inform us on the strength and direction of 
the component of the Galactic magnetic field perpendicular to the line-of-sight, they do not 
encode any information related to the component parallel to the line-of-sight. The latter can 
be provided by Faraday rotation studies (Sec. 2.1.2) or rotation measure synthesis (Sec. 2.1.3). 



2.1.2 Faraday Rotation of Compact Sources 



When a linearly polarized electromagnetic wave propagates through a region of free electrons 
permeated by a magnetic field (e.g., a magnetized plasma such as the ISM), its plane of 



polarization rotates. In the quasi-longitudinal (QL) approximation (Ratcliffe, 1959), which 
assumes the field direction is roughly along the line-of-sight, the Faraday rotation angle (in 
radians) at wavelength A (in meters) is ^' ~ RM, where RM, the "rotation measure," is 



RM = 0.812 j TicB-dl (radm"2) 



(1) 



with He the electron density in units of cm~'^, B the magnetic field in units of fiG, and dl the 
incremental pathlength from the source to the receiver in units of pc. At radio frequencies, with 
typical Galactic magnetic field strengths and electron densities, the QL approximation can be 
shown to hold for fields which are essentially perpendicular to the line-of-sight. Consequently, 
Eq. ([1]) is actually valid for most ISM applications. 

If we assume that at all wavelengths the polarized emission from a given source is emitted 
at the same initial polarization angle ipo, and that the radiation is only affected by Faraday 
rotation, the detected polarization angle is given by the linear relationship ip = (fo + X'^ RM. 
Measurements of ip at multiple wavelengths can therefore determine the RM for a given source 
as the slope of the graph of ip vs. A^. In practice, at least three frequency channels are needed 
to properly do this exercise. If many more are available, RMs may also be calculated by using 



a RM transfer function (RMTF) (see Brentjens fc de Bruyn 2005 and Sec. 2.1.3) 
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The simplicity with which RMs can be determined, coupled with the significance of their 
signs, makes them powerful tools with which one can probe the ISM magnetic field. Pulsars 
and extragalactic sources (EGS) emit linearly polarized radiation and are typically used to 
study the Galactic magnetic field. Since it is possible to estimate their distances to us, we can 
work backwards and find out what the field looks like along lines-of-sight to them, provided 
we make some assumptions about the electron density distribution within the Galaxy. 

When pulsars are used, studying the dispersion of the emitted pulses provides information 
to replace some of these assumptions. Different frequencies propagate at different group veloc- 
ities in the interstellar medium. As a result, a pulse emitted by a given pulsar is not received 
at the same time at different frequencies. The arrival time r at frequency v is such that 



dr 
d^ 



1 



oc 



rie dl , 



(2) 



which is a measurable quantity. Since the positive proportionality constant is a combination of 
known physical constants, it is possible to compute the value of the integral, called dispersion 
measure (DM), and therefore to estimate the integrated electron density along the line-of-sight. 

By measuring rotation measures for the highest density of sources possible, and over a 
wide range of Galactic longitudes, one can thus accurately reconstruct the large-scale Galactic 
magnetic field in the disk. Much effort has been put into these studies, and experiments such 



as the Canadian Galactic Plane Survey (Taylor et al. 2003) and the Southern Galactic Plane 



Survey ( McClure-Griffiths et al. 2005 Haverkorn et al. 20061 have now observed hundreds 



of square degrees near the Galactic plane with arcminute resolution. 



2.1.3 Rotation Measure Synthesis 

Rotation measure synthesis, also called Faraday tomography, is a novel technique for 3-D 



mapping of RM structures in the ISM. The idea was described by Burn ( 1966 1 decades ago, but 



only now are spectropolarimetric capabilities advanced enough that we can use this technique 



(Brentjens &: de Bruyn 2005). The Faraday depth of a rotating medium is defined as 



= 0.812 / neB • dl (radm" 



(3) 



where the integral runs over a particular segment of the line-of-sight, from x\ to X2, and rig, B, 
and dl have the same meaning and units as in Sec. 2.1.2[ In general, ^ RM. In addition, one 
can define the observed polarized flux density function as -P(A^) = W{}?) P(A^), where -P(A^) 
is the complex polarized surface brightness of the source as a function of wavelength, and 
W{}?) is a sampling function which is non-zero only at the wavelengths where measurements 
are taken. Fourier transformation of the polarized surface brightness as a function of gives 
F{(j)), the complex polarized surface brightness per unit Faraday depth (j) through 



+00 



P(A2 



(4) 



Thus, provided one can account for the effect of the sampling function, which is possible 



when many wavelengths are observed (Brentjens Sz de Bruyn, 20051, and make reasonable 



assumptions as to the definition of P for A < (Burn 19661, RM synthesis provides a 3-D 
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map of polarized intensity, where the third dimension is Faraday depth. However, Faraday 
depth does not necessarily increase monotonically with distance, as significant changes in field 
direction, electron density, or magnetic field strength will also impact its value, and there is 
therefore no direct distance information in these maps. Morphological correspondences with 
structures indicating a dramatic change in the magneto-ionic medium along the observed line- 
of-sight, such as spiral arms, supernova remnants, or external galaxies, can provide this missing 
link. Since the traditional way of measuring RMs only provides one RM value for each line- 
of-sight, and none where depolarization is severe, this new technique provides a tremendous 
increase in the amount of obtainable information from radio spectropolarimetric surveys. 



2.1.4 Other Techniques 



Polarized Dust Emission Non-spherical dust grains can partially align under the infiuence 
of an external magnetic field, which causes the emission from interstellar dust at far-infrared 
and submm wavelengths to be polarized. Although tremendous progress has been made in 



recent years (see, e.g., the review by Lazarian 20071, the grain alignment mechanisms are not 
yet fully understood. As the degree of polarization of this emission depends not only on the 
magnetic field strength but also on dust properties and alignment mechanisms, observations 
of polarized emission from interstellar dust do not currently allow the magnetic field strength 
to be directly inferred. However, the magnetic field direction averaged over the line-of-sight 
can be measured, as the polarization angle is perpendicular to the component of the magnetic 
field in the plane of the sky. This method has so far been mostly used in denser regions 



of the Galaxy, such as molecular clouds and cores (see, e.g., Vaillancourt , 2007 1 . In these 



environments, a rough estimate of the strength of the magnetic field can be obtained with the 



Chandrasekhar-Fermi method ( Chandrasekhar Sz Fermi 1953) by comparing the dispersion 



in polarization angle to the turbulent velocity dispersion. Numerical simulations suggest that 



results from this technique should be corrected by factors of order unity (Heitsch et al. 2001 1. 



Starlight Polarization When non-spherical dust grains are (partially) aligned with the Galac- 
tic magnetic field, the dusty region they are part of acts like a polarizer. As absorption of 
radiation by dust grains is more efficient in the direction of their long axes, initially unpolar- 
ized starlight illuminating the region becomes polarized in the direction of the grains' short 
axes, i.e., mostly parallel to the Galactic magnetic field (see Lazarian, 2007 for a more detailed 



discussion). By observing the polarization of the remaining starlight, inferences can therefore 
be made about the Galactic magnetic field required to produce the observed polarization close 



to the Sun, or, more generally, in regions of low dust extinction (see, e.g., Fosalba et al. 2002 1 



Zeeman Splitting In the presence of a magnetic field, degeneracies in atomic or molecular 
energy levels can be lifted. The energy separation between the previously degenerate Zeeman 
sub-levels is proportional to the strength of the magnetic field. As a result, a sufficiently 
large field must be present for the corresponding separation in a given radio-frequency line 
for a given species of atom or molecule to be detectable. In general, the magnitude of the 
interstellar magnetic field (~ 10 /iG) is not strong enough to generate detectable Zeeman 
splitting (see, e.g., Davies et al. 1960[ ) outside of localized interstellar clouds where the field 
is at least ~ 50 /iG (see, among others, the review by Davies, 2007 1 . 
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That being said, recent studies of Zeeman splitting from OH masers in massive star forming 



regions (Fish et al. 2003) and within molecular clouds (Han Sz Zhang 20071 have suggested 



that the fields within these regions are influenced by, and even reminiscent of, the large-scale 
Galactic field. While these results seem to be consistent with observations from other methods, 
more data is required before it can be considered a new way of constraining the Galactic field. 

Although Zeeman splitting can provide information on the strength of the Galactic mag- 
netic field, it does not encode information related to its direction. Observations of the Gold- 



reich-Kylafis effect (Goldreich Sz Kylafis 19811, which leads to a small linear polarization of 



interstellar radio-frequency lines, could complement Zeeman constraints from molecular clouds 
and star forming regions by providing the direction of the sky-projected component of the field. 

UHECR Multiplets The deflection angle of a ultrahigh energy cosmic ray (UHECR) is inde- 
pendent of the Galactic electron density, and is proportional to its charge Z and the transverse 
magnetic field, while being inversely proportional to its energy E. Thus, UHECRs originating 
from the same source will have their arrival directions appear, roughly, as a string on the 
sky, with the cosmic rays with higher value of Z/ E further from the true source direction. 
By detecting these multiplets, one can therefore estimate the direction of the source and the 



strength of the integrated transverse magnetic field along the line-of-sight (Jansson et al. 



2007 1 . This offers a direct way of constraining models of the Galactic magnetic field without 



having to rely on a model for the Galactic electron density (on this topic, see also Sec. 3.3). 



Although using this technique will have to wait until a bigger UHECR data set is available. 



it is worth noting that the Pierre Auger collaboration (Abraham et al. 2007), after 1.2 years of 



observations, found correlations between 28 cosmic ray events above 60 EeV and the positions 
of nearby AGNs. In ten years, over 200 events at that energy, including many multiplets, can 
therefore be expected in the Auger data, which could lead to potentially interesting constraints. 



2.2 Current Knowledge 

For practical reasons, the Galactic magnetic field is usually discussed as consisting of two 
components: the large-scale and the small-scale field. The large-scale, or "regular," magnetic 
field is predominantly azimuthal in the plane of the Milky Way, as clearly shown by the 



direction of starlight polarization by interstellar dust (Heiles 2000) and confirmed with the 



lowest-frequency maps in the WMAP 3- year polarization data (Page et al. 2007). The small- 



scale, or "random," component is of comparable strength as the regular field, exhibits structure 



over a wide range of physical scales, and is presumably turbulent (see discussion in Sec. 2.2.2 ). 



The average strength of the total magnetic field, i.e., regular and random components 
together, can be derived from observations of the Galactic synchrotron emission, as discussed 
in Sec. 2.1.1 Under the assumption of equipartition, the total magnetic field strength is 



estimated at 6 it 2 /xG in the solar neighborhood, and is found to increase to 10 it 3 /^G 



at a Galactocentric radius of 3 kpc (Beck, 2001). The equipartition assumption used in 



these estimates has been validated by measurements of the local 7-ray density, which yields 
comparable field strengths as a function of radius (Strong et al. , 2000). These results are 



consistent with those of Heiles (1996), who combined various methods to conclude that the 



total magnetic field strength is about 6 /uG in the spiral arms and 4 ^uG elsewhere. Near the 
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Galactic center, the magnetic field strength increases to several mG ( Yusef-Zadeh et al. , 1996 1 



In Sec. 2.2.1 we review the properties of the regular component of the Galactic field, 



whereas our knowledge of the random field is presented in Sec. |2.2.2 A Galactic halo field is 



also required to explain observations at low radio frequencies, which we discuss in Sec. 2.2.3 



2.2.1 The Large-Scale Magnetic Field in the Galactic Disk 

The local regular field strength estimated from pulsar RMs combined with their dispersion 



measures (DMs) is about 1.4 /uG (see, e.g.. Rand Sz Lyne, 1994), while estimates from radio 
synchrotron intensity give higher values of ~ 4//G. Because of uncertainties and assumptions 



in both methods (Beck et al. 20031, these results are not necessarily in contradiction. 



The large-scale magnetic field is observed to be concentrated in the Galactic disk (see. 



e.g., Simard-Normandin Sz Kronberg 1980; Han Sz Qiao 1994), and is generally assumed 



to be closely aligned with the Galactic spiral arms, consistent with observations of external 



galaxies (Beck, 2001 ). One feature that is critical in the determination of the dynamo mode(s) 
operating in the Galaxy is the presence of magnetic field reversals, regions of magnetic shear 
where the field direction changes by roughly 180°. One large-scale field reversal unquestionably 
exists between the local arm and the Sagittarius-Carina spiral arm of our galaxy, in the 



first quadrant (Simard-Normandin Sz Kronberg 1979 Thomson &; Nelson 1980). There are 



indications that this reversal slices through the Sagittarius-Carina arm near i = (Weisberg 



et al. 2004 Brown et al. 2007), suggesting a significantly different inclination of the field 



from the spiral arms, at least in this region. The existence of additional reversals is highly 
controversial (see, e.g. , |Vallee| |2002[ [Weisberg et"al| |2004! 'Han et al.| |2006| [Brown et al^ [2007[) , 



and recent modeling based on KM data (Noutsos et al. |2008, Men et al. , 2008 ), as well as on a 



combination of Faraday rotation and diffuse Galactic synchrotron emission measurements (see 



Sec. 3.3), all conclude that the current data is not sufficient to conclusively distinguish between 



several large-scale magnetic field models predicting differing numbers and locations of reversals. 
It is worth pointing out that magnetic field reversals have yet to be observed in external spiral 



galaxies (Beck, 2007). It is uncertain if this is a consequence of observation techniques and 



resolution, or a true indication of some uniqueness in the properties of the Galaxy. 

Figure [T] illustrates the directions of the large-scale Galactic magnetic field superposed on 



the electron density model of Cordes Sz Lazio (2002). Both universally accepted (solid arrows) 
and still debated (dashed arrows) directions are shown (see caption of Fig. [T] for more details). 

2.2.2 The Small-Scale Magnetic Field in the Galactic Disk 

The ratio of the regular to the small-scale Galactic magnetic field strength is estimated from 
synchrotron emission to be approximately 0.6 ~ 1 (see, e.g., Phillipps et al. 1981 Heiles 1995 ), 



in agreement with estimates from the variance in pulsar RMs (Rand Sz Kulkarni, 1989). 



As the energy density in the magnetic field is comparable to the (turbulent) gas density 



(Boulares Sz Cox 1990), interstellar turbulence is expected to be magnetic, although details 



about its structure, strength and variability are still largely unknown. Interstellar turbulence 
is traditionally measured through power spectra of either velocity fluctuations in emission or 
absorption lines, or density variations. Similarly, one can get an idea of the role of the magnetic 
field by studying power spectra of RM measurements. To avoid the problem of depolarization 
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along the line-of-sight, this is usually done using RMs of pulsars and polarized extr agalactia 
sources. The corresponding data points being irregularly spaced, the computation of a power 
spectrum is usually replaced by a structure function analysis, with the second order structure 
function of RM given by DRM(dx) = (RM(x+dx) —RM{x))x, where dx is the distance between 
two sources, and denotes averaging over all positions x. These analyses show that RM 
fluctuations exist over a wide range of scales and exhibit a power law behavior (see, e.g., 
Minter &; Spangler 1996). The largest observed fluctuations are on scales of 50 - 100 pc when 



averaged over the sky (Rand &: Kulkarni, 1989 Ohno &: Shibata 19931, and of order 1 pc in 
the spiral arms ( Haverkorn et al. 2008 1 . At least part of this structure is caused by magnetic 
field fluctuations as opposed to variations in the electron density (Minter & Spangler, 1996 1 . 




Figure 1: Directions of the large-scale Galactic magnetic field as viewed from the North Galactic 
pole. The grey scale represents the electron density model of jCordes fc Lazio (20021, Ql-4 
indicate the four Galactic quadrants, and the asterisk shows the location of the Sun. Solid 
arrows correspond to universally accepted field directions, while single-ended dashed arrows 
show field directions as supported by the Canadian or the Southern Galactic Plane Survey 



(CGPS and SGPS, respectively), though not necessarily universally accepted (Brown et al. 



2003b[p007| . Finally, double-ended arrows indicate regions remaining highly debated with 



either no data or no clear indications currently available from the CGPS or the SGPS. 



2.2.3 The Magnetic Field in the Galactic Halo 

Significant synchrotron emission is observed at high Galactic latitudes, indicating the presence 
of a magnetic field in the Galactic thick disk or halo. The regular magnetic field has a scale 
height of about 1.5 kpc and is an order of magnitude weaker than the disk field, as derived 
from pulsar RMs (Han & Qiao 1994 1 and supported by extra-Galactic RMs (Brown et al., 
in preparation). In the thin disk of the Galaxy, the field points in the same direction above 
and below the Galactic plane (Frick et al. 2001 1 . However, there is evidence for large-scale 
reversals of the field's direction in the thick disk/halo, which are likely caused either by a 
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dipole structure of the halo magnetic field (Andreasyan &; Makarov 1988 Han et al. 19971, 
or by local structures that appear large on the plane of the sky. 



3 Global Constraints on the Magnetized ISM 



As reviewed in Sec. |2.1[ numerous techniques can be used to probe the structure of the 
Galactic magnetic field on various scales. However, most of these studies do not lead to direct 
constraints on the latter, but on the magnetized interstellar medium, of which the Galactic 
magnetic field is only one component along with, e.g., thermal and cosmic ray electrons. Given 
a model for the physical properties of these other components, each technique can then be 
used to characterize the field in our galaxy. These models are usually themselves derived from 
other observations, and their reliability is often limited, which, in turn, weakens the constraints 
one can infer for the Galactic magnetic field. An example can be found in the controversy 



surrounding the existence of field reversals in the outer parts of the Galaxy (see Sec. 2.2.1). 

Although this approach has several advantages, in particular insofar as it offers the possi- 
bility of comparing independent constraints on the various components of the field (see, e.g.. 



Crutcher et al. 20031, it fails to account for the possible degeneracies between properties of 



the field and characteristics of other components of the magnetized interstellar medium. This 
can be avoided by simultaneously fitting multiple observations, which also has the advantage 
of leading to a unique model able to describe a diverse set of data. Such a model is all the 
more desirable since the same physical quantities are involved in various probes of the field. 

There are multiple data sets one can combine to build such a model, but one first needs a 
simulation code able to generate the corresponding theoretical maps from a set of parameters 
describing the relevant components of the magnetized interstellar medium. The Hammurabi 



software (Waelkens et al. 2008) was developed to provide such a tool for the Galactic syn- 



chrotron emission (total and polarized), free- free emission, and the effect of Galaxy-induced 
Faraday rotation on the emission from extra-galactic sources. In this section, we report on 



the efforts by two groups: Jansson et al. (2007) and Sun et al. (2008). The former used Ham- 



murabi to constrain the magnetized interstellar medium with the three-year WMAP 23 GHz 



polarization data (Page et al. , 2007). Sun et al. (2008) then performed an extensive analysis 



taking into account the total synchrotron emission observed by Haslam et al. ( 1981 1982 ), the 



1.4 GHz synchrotron map of WoUeben et al. (2006), the WMAP three- year free- free emission 



template (Hinshaw et al. 2007), and the rotation measures of Brown et al. (2003b 2007) 



3.1 ISM Modeling 

Assumptions about three components of the magnetized ISM, namely, thermal electrons 
(Sec. 3.1.1), cosmic ray electrons (Sec. 3.1.2), and the Galactic magnetic field (Sec. 3.1.3), 
are necessary to make it possible for Hammurabi to generate the simulated maps mentioned 
in Sec. |3j Comparison to data then enables Jansson et al. (2007) and Sun et al. (2008) to con- 
strain the parameters of these models by making sure that they can reproduce their selected 
sets of observations. We review each of these models and their associated parameters below. 



Some of the constraints obtained by the Jansson and Sun groups are reviewed in Sec. 3.3 
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3.1.1 Thermal Electron Density and Temperature 



Hammurabi uses the NE2001 model developed by iCordes k Laziol (120021 120031) to describe 



the Galactic thermal electron density rig. This model, based on dispersion measures, contains 
both a thick disk and a thin disk, as well as spiral arms and some local ISM features. 



Sun et al. (20081 combined the results of Reynolds et al. (19991, Peterson k Webber (20021 



and Quireza et al. ( 2006 1 to show that the temperature distribution of the thermal electrons 
is well described by Te{r)/K = 5780 + 287 r — 526 \z\ + 1770 z^, where r and z (cylindrical 
coordinates) are in units of kpc. Moreover, the free-free emission intensity is proportional 
to (jT-e^), where the average is along a given line-of-sight, which is a priori different from 
the quantity (ng)^ that can be inferred from the NE2001 model. This discrepancy can be 
accounted for by the filling factor fe{z) 



(ue) /{ue ) provided by Berkhuijsen et al. (20061 



3.1.2 Cosmic Ray Electron Density 



Models of the cosmic ray electron density such as those used by Page et al. (2007) and Sun 



et al. (2008) are composed of a simple spatial distribution characterized by a scale height and 



a scale length, and also assume a power law energy distribution for the electrons. The latter 
has so far been sufficient to model foregrounds in the context of CMB data analysis [see, e.g.. 



Tegmark k Efstathiou (1996), de Oliveira-Costa et al. (2008), Dunkley et al. (2008b), and 



Bottino et al. (2008)]. However, simulation codes, such as GalProp, lead to models much 



closer to the state-of-the-art in cosmic ray modeling (as in Strong et al. 2007) and might have 



to be used in conjunction with more sophisticated analyses in the future (we refer the reader 



to the companion paper by Dunkley et al. (2008a) for details on this issue 



The Sun et al.| ( |2008[ ) model assumes a scale height (respectively, a radial scale length) of 
1 kpc (respectively, 8 kpc), whereas the normalization is chosen to reproduce observations at 
r = r0 (see, e.g.. Strong et al. 2007). Furthermore, a truncation is introduced for |z| > 1 kpc 



to explain the low synchrotron emission observed at high Galactic latitudes. Finally, the 
energy distribution is given by the power law N{E) oc E~p, where p is 3 for u > 408 MHz, 
and 2 at lower frequencies. No local features of the Galactic interstellar medium are modeled. 



3.1.3 Galactic Magnetic Field 



Any large-scale magnetic field model can be tested with Hammurabi. Jansson et al. (2007) 



looked at axisymmetric and bisymmetric spiral models, with symmetry or antisymmetry in 



the z direction, while Sun et al. (2008), after demonstrating the need to improve models found 



in the literature, considered an axisymmetric spiral field with a ring, an axisymmetric spiral 
field with a field reversed arm, and a bisymmetric spiral field. In addition, a halo magnetic 
field, necessary to fit high-latitude rotation measure profiles, and a random field, required by 



the synchrotron and RM observations, are included in the Sun et al. (2008) analysis. 



As for the cosmic ray electron model, no attempt has been made at modeling the Galactic 
magnetic field to explain local features such as the North Polar Spur. Although they occupy 
a large fraction of the sky, and are fascinating objects to explore on their own (see, e.g.. 



the Loop I model developed by Wolleben 2007), they are beyond the scope of a first global 



description of the diffuse Galactic emission, and are therefore not included in Hammurabi. 
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3.2 Simulation Outputs 



Based on the model described in Sec. |3.1[ Hammurabi can generate Galactic simulated maps 
of synchrotron emission (total and polarized), free- free emission, and rotation measure (Fig.|2|. 
Structures like the North Polar Spur, or any of the other Loops, are not reproduced by these 
simulations since the magnetic field and cosmic ray models only reproduce the large-scale 
features of the Galaxy. However, local features in the electron density distribution, which are 



part of the NE2001 model (see Sec. 3.1.1), can be recognized in the rotation measure map. 



total mtensit}' {408 MHz; 



,ed intensity (22 GHz) 




7* < . 



I 0,12 log (mK) 



polarization angl 



rotation measure (KM) 





Figure 2: Mock observations of the 



Sun et al. 



( 2008 1 model. The top row shows full-sky maps of the 



total (left) and polarized (right) synchrotron emission, whereas the bottom row provides 
a map of the corresponding polarization angle (left) and a rotation measure map (right). 



Unlike in Sun et al. (2008 1, the small-scale magnetic field is modeled as a Gaussian random 



field. The simulated map of free-free emission is not shown here. 



3.3 Comparison to Observations 



Using a % minimization technique, Jansson et al. (20071 compare the polarized synchrotron 
maps produced by HAMMURABI to the three-year WMAP K-band (23 GHz) map smoothed 
to 4°. They find that the minimized is close to identical for all models of the Galactic 
magnetic field they consider (see Sec. 3.1.3). However, for each of them, large regions of 
the parameter space can be excluded, even though there are no obvious best-fit parameters 
as a number of combinations yield a value of very close to the minimum. In particular, 
the vertical and radial scale-lengths can take unphysically large values, which is explained by 
the drop in synchrotron emission due to the reduction in the relativistic electron density at 
large radial and vertical distances. Thus, knowing the polarized synchrotron emission is not 
sufficient to strongly constrain the parameters of these Galactic magnetic field models. 
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By combining synchrotron emission data (Page et al. 20071 with rotation measure infor- 



mation (Brown et al. , 2003b 2007) and knowledge of the full-sky free- free emission (Hinshaw 



et al. 2007), Sun et al. (2008) showed that this degeneracy can be at least partially lifted. They 



indeed demonstrated that increasing the vertical scale-length of the NE2001 thermal electron 
density, removing the unphysical cut in the cosmic ray electron density (see Sec. 3.1.2), and 
decreasing the strength of the halo magnetic field, leads to a better global fit to the observa- 
tions. The need for a thermal electron scale height about twice as large as that of the NE2001 
model was confirmed by an independent data analysis by Gaensler et al. (2008). 



Sun et al. (2008) were also able to find statistical evidence that the Galactic magnetic field 



is best fit by an axisymmetric spiral field with a ring located in the inner part of the Galaxy, 
thereby adding to the case that combining observations probing a broad range of physical 
processes is a powerful way of constraining the various components of the magnetized ISM. 



3.4 Future Developments 



In addition to the probes combined by Sun et al. (2008), other observations could already 
be used to improve our knowledge of the magnetized interstellar medium. An example is 
pulsar RM - DM plots (see, e 



Brown et al. 


2003b 


and Sec. 


2.2.1 



dimensional structure of the Galactic magnetic field. Hammurabi's design allows it to perform 
the corresponding mock observations, and the Square Kilometer Array (SKA) is expected to 



provide thousands of data points for these plots by 2020 (Beck Sz Gaenslerl 2004). 



The Sun et al. (2008) model does not contain any dust related inputs, and, as a result, 
cannot make starlight polarization or dust emission predictions. Dust emission mechanisms 
are difficult to model given the sophisticated and sometimes unknown physics associated with 
them. Moreover, models developed to match all available observations can lead to significantly 
different predictions for what the polarized emission from interstellar dust might look like (see 



Sec. 5.3). Adding a dust component to Hammurabi will therefore be a challenging task. 

Upcoming radio telescopes, such as the SKA, the LOw Frequency ARray (LOFAR), and the 
Atacama Large Millimeter Array (ALMA), associated to upcoming CMB experiments, such 
as Planck, will greatly increase the number and the quality of available observations. A tool 
like Hammurabi will allow us to make the most of these improvements, and, in particular, 
to generate valuable information regarding the structure and properties of the magnetized 
interstellar medium. This will be a welcomed source of inputs for the CMBPol data analysis. 



4 The Galactic Magnetic Field in the CMBPol Era 

4.1 Radio Surveys in the Next Decade 

Narrow frequency and spatial coverage have so far been the main limiting factors in studying 
the low frequency polarized emission from the Galaxy. This issue is expected to be addressed 
by large radio polarimetric surveys currently planned or underway, which will provide a wealth 
of new data by the time CMBPol flies. The S-PASS polarimetric survey (PI: E. Carretti) will 
observe the entire Southern sky in a 256 MHz frequency band centered on 2.3 GHz by 2009. 
The GALFACTS polarized continuum survey (PI: A. R. Taylor) will map the entire sky vis- 
ible from Arecibo in a 300 MHz frequency band centered on 1.4 GHz by 2012, whereas the 
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STAPS survey (PI: M. Haverkorn) will cover 500 MHz around the same frequency in the 
Southern sky. Radio polarimetric observations are also planned with the Australian SKA 



Pathfinder (Johnston et al. , 20071, which should be fully operational by 2013. Finally, the 



Galactic Magneto-Ionic Medium Survey (GMIMS, PI: M. WoUeben), a major project consist- 
ing of six independent radio polarimetric surveys, is collecting high-resolution (~ 30') all-sky 
polarization data in a near-continuous frequency band ranging from 300 MHz to 1.8 GHz. 



These and similar surveys will provide "RM grids" (Beck &; Gaensler 2004 1 gathering 



numerous measurements of the Galactic Faraday rotation of polarized emission from extra- 
galactic sources across the whole sky. Although EGSs have a small intrinsic RM component, it 
can be statistically eliminated by averaging the RMs of enough point sources in a given patch 
of the sky. The current average RM source density is about 1 RM source per 10 square degrees 
at high latitudes, and approximately 1 source per square degree in parts of the Galactic plane 



(Brown et al. 2003a 2007 1 . As a result, inferred Galactic RM maps are very low resolution 



over most of the sky. However, the radio surveys described above are expected to fill in the 
Galactic plane to the level of 1 source per square degree, as well as to increase the source 
density at high latitudes by a factor of 100. This will allow averaging of a sufficient number 
of RM sources to remove intrinsic RM contributions on scales of about 1°. Detailed modeling 
of the pitch angle of the Galactic magnetic field, the number and locations of large-scale field 
reversals, and the vertical component of the magnetic field will therefore soon become possible. 

Moreover, RM grids will enable detailed study of the turbulent field component over a wide 
range of scales and large parts of the sky. Given their expected high source density, it will 
become possible to construct a structure function every few square degrees on the sky, which 
will result in unprecedented detailed knowledge of the turbulent behavior of the magnetized 
ISM on degree scales. In particular, possible changes in the turbulence characteristics from 
the Galactic plane, via the disk-halo connection to the Galactic halo, should become visible. 

Finally, RM synthesis of polarized emission from diffuse synchrotron emitting gas in the 
Galaxy can be used to obtain a 3-dimensional map of Galactic RM structure. All of the 
above mentioned surveys have the capability to perform RM synthesis on their data, with 
various limitations, which will be used to improve modeling of the large-scale components of 
the Galactic magnetic field. As the technique allows separation of magnetic foregrounds and 
backgrounds, it will also be possible to study magnetic field structures in Galactic objects such 
as supernova remnants, HII regions, and planetary nebulae, as well as in external galaxies. 



4.2 The Contribution of Upcoming CMB Experiments 

To study the global character of the Galactic magnetic field through synchrotron emission 
and (de)polarization, observations over a significant part of the sky are needed, at a frequency 
high enough that Faraday rotation is negligible. WMAP provides polarization maps that 



allow study of the Galactic magnetic field through synchrotron radiation ( Page et al. 2007 1 



However, the polarized emission in the lowest frequency WMAP maps has a signal to noise 
ratio below 3 in about 55% of the sky when smoothed to 2°. The only polarized structures 
that are clearly visible in the WMAP maps are the Galactic plane, the region of high and 



regular polarization called the Fan region (Berkhuijsen et al. 1964 Brouw &; Spoelstra 19761 



and the North Polar Spur (Hanbury Brown et al. 19601, generally believed to be an old 
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nearby supernova remnant. Therefore, to be able to study the behavior of the synchrotron 
radiation away from the Galactic plane in the "typical" interstellar medium, higher sensitivity 
observations are needed. WMAP polarization maps show that although in both the Galactic 
plane and discrete structures the polarized intensity can go up to ~ 100 ^K, at intermediate 
and high latitudes the polarized intensity is primarily below 20 ^K. This will allow a nominal 
detection with Planck or nine years of WMAP observations, but with a typical signal-to-noise 
lower than 2, which will lead to an error in polarization angle determination of at least 1 rad. 
CMBPol, which could be 10 times more sensitive than Planck at the same frequencies, would 
have much higher signal-to-noise across the whole sky, but in particular in these intermediate 
and high Galactic latitude regions. This data would complement the information provided by 



the radio surveys described in Sec. 4.1, in particular with respect to the vertical component 
of the field, as well as the strength and structure of the halo field, both of which are necessary 
knowledge to test dynamo theories of the origin and evolution of Galactic magnetism. 

In addition, the resolution of CMBPol might be higher than that of Planck, which would 
significantly decrease depolarization. At the frequencies currently considered for CMBPol, 
there is no significant Faraday depolarization, but small-scale structure in the magnetic field 
will cause small-scale structure in polarization angles of the synchrotron emission, resulting 
in depolarization of the radiation within one telescope beam. If the satellite's resolution is 
not high enough, the polarized emission will be below its detection threshold, destroying all 
diagnostic capability of the data for Galactic magnetic fields. However, at high latitudes, 
turbulent cells are expected to be at ~ 100 pc scales (Dumke et al. 19951, which, at the dis- 



tances of interest, corresponds to a much larger angle than the resolution currently considered 
for CMBPol. As a result, depolarization is not expected to be severe there. On the other 
hand, the typical structure of the magneto-ionized medium in spiral arms could be on scales 
as small as 1 pc (iHaverkorn et al. 20081, which makes depolarization resolution dependent. 



In particular, higher resolution than Planck^s 33' beam at 30 GHz will be necessary to observe 
turbulent structures in the magneto-ionized medium in and close to the Galactic plane. 

The same holds for discrete polarized structures, such as supernova remnants (SNRs). 



Statistical studies of supernova rates (see, e.g., Li et al. , 1991) indicate that there should be 
many more SNRs in the Galaxy than are currently known. This "missing SNRs" problem 
is likely due to a selection effect of many radio surveys against large low-brightness SNRs 
as well as small remnants. However, low-brightness remnants might be observable through 
their magnetic field structure, which remains visible through polarized radio emission even 



when synchrotron intensities are too low to be observed ( Haver korn 2005). A sensitive all-sky 



survey like CMBPol with the same resolution as necessary to study the turbulent structure of 
the field in spiral arms therefore opens the possibility of detecting the magnetic field structures 
of old and faint SNRs, which could help solve the long-standing "missing SNRs" problem. 



5 Polarized Dust Emission (PDE) 

Over the course of the last 60 years, the linear polarization of starlight has been observed at 



wavelengths ranging from the ultraviolet (Clayton et al. 1992 Anderson et al. 1996 Martin 



et al. 19991 to the far-infrared (see, e.g., the review by Whittet 2003) along numerous lines- 



of-sight. The Archeops experiment first released detailed studies of the emission properties of 
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dust at CMB wavelengths based on observations of 20% of the sky at 353 GHz ( Benoit et al. 



2004 Ponthieu et al. 20051, and warned that the polarized dust emission would be a major 



contaminant of the primordial E- and B-mode signals unless removed properly. These results 



were largely confirmed by the WMAP mission (Kogut et al. 20071, although the fact that the 



highest WMAP frequency (94 GHz) is below 100 GHz made it easier to handle this emission 



in the polarization analysis (Page et al. 2007 Gold et al. 20081. 



The Planck satellite will soon be on its way to L2, where it will map the sky at frequencies 
ranging from 30 GHz to 857 GHz, with 5 of its bands covering frequencies higher than 100 GHz. 
Compared to WMAP, Planck will have twice as many frequency channels, 3 times higher 
resolution, and 10 times higher sensitivity, all of which will, by design, be of tremendous help 
to distinguish between polarized foregrounds and candidate inflationary signals. 

However, a number of assumptions had to be made to reach this "optimal" configuration. 
In particular, the dust emission, which will be brighter in the 5 Planck frequency bands 
mentioned above than in any WMAP map, was assumed to have the same polarization degree 
at all wavelengths, with the 857 GHz observations providing a way of determining, among 
other information, this polarized fraction. Although so far consistent with the degrees of 



polarization inferred by the Archeops and WMAP groups (Benoit et al. 2004 Kogut et al. 



20071 at, respectively, 353 GHz and 94 GHz, there is no strong observational evidence that 



this will hold at all Planck frequencies. If the level of polarized dust emission varies as a 
function of frequency, it will reduce our ability to remove this emission from the Planck maps. 
A CMBPol mission is likely to have between 4 and 6 frequency bands in which the brightest 
signal will be dust emission, and might therefore have to deal with the same problem. 



In this section, we follow Draine Sz Fraisse (2008), and show how our current knowledge of 



the extinction and polarization of starlight can be used to predict the degree of polarized dust 
emission one could observe at CMB wavelengths with a high-resolution and high-sensitivity 
mission such as Planck or CMBPol, as well as its potential dependence on frequency. 



5.1 Observed Extinction and Polarization of Starlight 

Although observations of the extinction and polarization of starlight exist at numerous wave- 
lengths [see references in the first paragraph of Sec. [5} as well as Martin et al. (1992), and, e.g.. 



the review by Fitzpatrick (2004)], they are well fit by simple parametric functions. Fitting the 



average Galactic extinction curve is discussed in Fitzpatrick (1999) for 0.1 < A/^m < 2.65 



whereas Anderson et al. (1996) showed that the fractional polarization curve p{X) is well 
described by a Serkowski law (Serkowski 1973) for 0.15 < X/fim < 1.39. Finally, for 
1.39 < X/fim < 5, p{X) oc A~^'^ is a good match to the observations by Martin et al. (1992). 



It is therefore possible to represent the observed extinction and polarization of starlight by 
smooth continuous functions in the range of wavelengths for which data is available. 

The fact that we observe starlight polarization implies that dust grains cannot be spherical, 
and that they have to be at least partially aligned with the Galactic magnetic field. However, 
the data currently in hand does not provide direct constraints on the shape of grains, their 
size distribution (s), or their alignment function(s). As a result, modeling interstellar dust 
requires relying on theory to predict the general behavior expected for the latter. One can 
then assume a shape and try to fit the observed extinction and polarization of starlight. Of 
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course, this process does not guarantee that there is a unique good fit as only a few shapes can 
be tested in a reasonable amount of time. But the hope is that, given the many constraints 
imposed on candidate models, if several of them can accurately reproduce observations of the 
extinction and polarization of starlight over a wide range of wavelengths, they will lead to 
similar predictions for what the polarized dust emission should look like at CMB wavelengths, 
at least as long as the non-spherical interstellar grains are made of the same material(s). 

5.2 Modeling Interstellar Dust 

Dust grains are believed to be made of the numerous elements that have been depleted from the 



gas phase (Jenkins 2004 1 . However, observed spectroscopic features in the Galactic extinction 



curve and in observations of emission by interstellar dust in the mid-infrared indicate that the 
bulk of the mass in interstellar grains is in the form of amorphous silicates and carbonaceous 
materials, among which a population of Polycyclic Aromatic Hydrocarbons (PAHs) is required 



to explain several emission features (see, e.g., the review by Draine, 20031. 



In order to explain the polarization observed in some of the spectroscopic features that have 
been studied in the far-infrared, amorphous silicates have to be present in non-spherical dust 



grains aligned with the Galactic magnetic field, B (see, e.g., Whittet 20031. As of today, there 



is however no compelling evidence that the same is true for carbonaceous materials (Chiar 



et al. 2006 ) . Finally, PAHs are expected to have their angular momenta only slightly aligned 



with B, and a PAH's principal axis of largest moment of inertia will be only partially aligned 



with its angular momentum (Lazarian Sz Draine 20001. Furthermore, PAHs are thought to 
contribute only a small fraction of the total emission by dust at frequencies v < 10^ GHz. As 
a result, PAHs should contribute negligibly to the polarized emission at those frequencies. 



5.3 Far-Infrared Emission from Interstellar Dust 



Based on the considerations discussed in Sec. 5.2 Draine & Fraisse (2008) fit the observed 



extinction and polarization of starlight with dust models consisting of (1) a population of 
(nearly) randomly-oriented PAHs, (2) oblate spheroidal silicate grains, and (3) graphite grains 
assumed to be either spheres (in which case their emission is not polarized) or oblate spheroids 
with the same axial ratio b/a of their equatorial to their polar axes as the silicates. 



Table 2: Draine & Fraisse (20081 Dust Models 



Model number 


(b/a) silicate 


(b/ o) carbon 


1 


1.4 


1.0 


2 


1.4 


1.4 


3 


1.6 


1.0 


4 


1.6 


1.6 



Table [2] provides the values of this ratio for each component of the Draine Sz Fraisse (2008 ) 
dust models, while Fig. [3] shows the corresponding predicted far-infrared emission spectra, as 
well as its expected degree of polarization as a function of wavelength, under the assumption 
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that all dust grains are heated by the Mathis et al. (1983) Galactic radiation field, and that 
the observer is looking in a direction perpendicular to the Galactic magnetic field. 

Despite the wide range of axial ratios probed by the 4 models listed in Tab. |2] for both 
amorphous silicate and carbonaceous grains, Draine &: Fraisse (20081 showed that they all 
provide adequate fits (not shown here) to the observed extinction and polarization of starlight, 
and give far-infrared emissions compatible with the DIRBE measurements (see Fig. [3| . The 
emission predicted at CMB wavelengths is also remarkably similar between all models. 
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Figure 3: {Left) Far-infrared emission spectra and {Right) corresponding degrees of polarization cal- 
culated by Draine & Fraisse (20081 for the models in Tab. |2] On each plot, the vertical 
dotted lines indicate 4 possible CMBPol frequencies ranging from 100 GHz to 850 GHz. 
The symbols on the emission spectra represent the HI correlated emission as measured by 



DIRBE at high Galactic latitudes (Dwek et al. 1997 Arendt et al. 19981 



However, the expected degree of polarization of the far-infrared emission from interstellar 
dust and its wavelength dependence are strong functions of the class of models one considers. 
When both silicate and carbonaceous grains are non-spherical and partially aligned with the 
Galactic magnetic field, the fractional polarization of the dust emission is essentially indepen- 
dent of frequency at CMB wavelengths, and of order 10%. But if only the emission by silicates 
contributes to the polarized signal, the polarized fraction of the dust emission can vary by a 
factor of 1.4 between 850 and 100 GHz, from around 11% to over 15%. For a given class of 



models, the level of polarization calculated by Draine &: Fraisse (2008) is not very sensitive to 



the value of the ratio h/a considered, which makes their predictions relatively robust. 



5.4 Polarized Fraction of Dust Emission for CMB Forecasts 

The polarized fractions shown in Fig. |3]may seem high in light of the degrees of polarization 
reported by Archeops ( [Benoit et ah| |2004[ ) or WMAP ( |Kogut et ah] |2007[ ) for the diffuse 
Galactic emission. However, there should be very few lines-of-sight in the sky close to perpen- 
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dicular to the Galactic magnetic field, so the average numbers these groups report, ranging 
from 4% to 7%, should indeed be smaller. That said, on a given sightline, the ratio of the 
observed polarization to the polarization calculated for any of the models in Tab.|2]should not 
be a function of frequency. Therefore, if one believes that the Galactic magnetic field models 
used by Archeops and WMAP are accurate enough to get a good estimate of the polarized 
fraction of dust emission of 5% at 100 GHz, we can easily rescale the Draine &: Praisse (20081 



models to reproduce this observation. All of them are then consistent with a polarized fraction 
between 3% and 5% at all frequencies CMBPol would potentially observe. 

Depending on the precision to which the polarized fraction 11 of dust emission needs to be 
determined in order to achieve a mission's science goals, such as detecting B-modes, several 
scenarios might therefore have to be considered in a CMB forecasting pipeline: 

• n = x% and is independent of frequency for v < 900 GHz [as is approximately the case 



for models 2 and 4 of Draine &: Fraisse (20081], 

H is a decreasing function of frequency, smoothly varying from H = x% at 



100 GHz 



to ~ 0.7 x% at ~ 850 GHz [as in models 1 and 3 of Draine fc Fraisse (2008)], 



where x is chosen either as having a value corresponding to one of the models shown in 
Fig. |3] or from an estimate obtained by the Archeops or the WMAP groups. In the former 
case, one should keep in mind that the chosen value will likely be an overestimate of the 
observable degree of polarization, whereas in the latter case, the polarized fraction might be 
underestimated compared to its true value due to the coarse resolution and low sensitivity of 
the Archeops and WMAP instruments compared to the Planck or the CMBPol detectors. 



5.5 Constraining Dust Models with CMBPol 



With its many channels and its high resolution, high sensitivity detectors, CMBPol should 
be able to average the polarized fraction of dust emission on much smaller patches than 



WMAP (Kogut et al. 2007), with potentially higher precision, and to detect variations of the 



degree of polarization between frequency channels if it exists, and if the polarized fraction of 
dust emission is high enough. How well this can be achieved will of course depend on various 
factors, including the quality of the Galactic magnetic field models that we will be able to put 
together by the time the CMBPol data is analyzed. As a result, it is hard, if not impossible, 
to know what CMBPol will bring to the debate on what interstellar dust grains look like, and 
how they behave. However, we can envision two main game-changing scenarios: 



Consistent observations of low polarized fractions of dust emission (typically. 



< 



at 



all CMB frequencies would suggest that the Draine & Fraisse ( 2008| ) models in which 
only silicate grains contribute to the polarized emission are not an obvious good match 
to the measurements. In this event, we might have to conclude that carbonaceous grains 
are significantly non-spherical and somewhat aligned with the Galactic magnetic field 



(as in their models 2 and 4), a fact for which we have currently no evidence (see Chiar 



et al. (2006) and discussion in Sec. 5.2). Alternatively, it could mean that the adopted 



opacities of the silicate material have to be reduced (in a wavelength-dependent fashion) , 
while those of the carbonaceous grains should be correspondingly increased. 
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• However, if the polarized fraction of dust emission is found to be significantly decreasing 
when frequency increases, and potentially higher than what the Archeops and WMAP re- 
sults suggest, this would be another indication that amorphous silicate grains and car- 
bonaceous grains indeed form separate populations, with the polarized emission being 
dominated by the former. This would pose a challenge to theories of grain alignment. 

We also have to mention the possibility that a behavior completely different from the results 
shown in Fig. [s] could be observed. Although the models of Draine Sz Fraisse (2008 1 appear to 
be consistent with all current observations, our understanding of interstellar dust is far from 
perfect. Such observations would therefore likely lead us to revisit some of the assumptions of 
our models for interstellar dust, opening a new window on dust properties and dust physics. 



5.6 Tracing PDE with Starlight Polarization 



The Draine Sz Fraisse ( 2008 1 results predict the polarized fraction of far-infrared and submil- 



limeter dust emission, and its dependence with wavelength, for lines-of-sight perpendicular to 
the Galactic magnetic field. Going from observed degrees of polarization to a function directly 



useable to constrain dust models (see Sec. 5.5), and, in particular, theories of grain alignment 



(see Sec. 6.3), therefore requires as good a large-scale magnetic field template as possible. 



As mentioned in Sec. 2.1.4, optical or near-infrared starlight polarization gives us the sky- 
projected direction of the Galactic magnetic field integrated along lines-of-sight going through 
regions of low dust extinction (typically. Ay or Ak < 5 mag). As a result, not much can be 
learned from this technique when it comes to studying the large-scale magnetic field in the 
Galactic disk, since it cannot probe the field in or near the Plane outside of the solar neigh- 
borhood. However, at intermediate and high Galactic latitudes, |6| > 10°, dust extinction is 
relatively low (see, e.g., Burstein &; Heiles 1982 and references therein), with Ay — 0.03 mag 
near the Galactic poles (the estimated error on this number is of order the extinction itself). 
Since these latitudes correspond to the regions of main interest to the CMB community, useful 
magnetic field templates can therefore be constructed from starlight polarization observations. 



Heiles (2000) compiled a catalog, the largest to date, of the polarization of starlight from 
nearly 9300 stars. However, among those, Fosalba et al. (2002) only found ~ 1400 stars with 



reliable distance and extinction estimates lying sufficiently far from the Galactic plane, i.e., 
with |6| > 10°, to be useful for CMB purposes. Moreover, the vast majority of these stars, 
about 1300, are less than 1 kpc away from the Sun. Despite this lack of information outside 
of the plane, starlight polarization data has already been useful to CMB polarization studies. 



Page et al. ( 2007 ) and Dunkley et al. ( 2008b ) both used the Heiles ( 2000 ) catalog to inform their 



fit to the dust emission component of their signal, and Kogut et al. (2007) showed that using the 



polarization angle provided by starlight polarization observations leads to a (marginally) better 
fit to this component than a model based on WMAP^s low frequency polarized synchrotron 
observations. Since both starlight polarization and polarized dust emission are processes 
closely related to how grains align with the interstellar magnetic field and to the optical 
properties of interstellar dust grains, this improvement is not totally unexpected. 

The scarcity of optical and near-infrared starlight polarization data has led several groups 
to try and remedy the situation, especially in view of the accuracies now achievable (see, e.g.. 



Carciofi et al. 2007). An example is the Galactic Plane Infrared Polarization Survey (GPIPS), 
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which aims at measuring nearly half a million i7-band stellar polarizations for H down to 



12 mag across a 72 square degree region of the inner Northern Galactic plane by 2011 (Clemens 



et al. , 2007). But neither the GPIPS nor other surveys currently planed or underway (see, 



e.g., Magalhaes et al. 2005; Nishiyama et al. 2008) will lead to a significant increase in the 



number of available data points in the diffuse interstellar medium at high Galactic latitudes. 
Given the implications this new data would have for our understanding of dust processes, but 
also for our ability to separate polarized dust emission from a potential gravity wave signal, 
we believe that acquiring it should be a high priority goal achieved by the time CMBPol flies. 



6 Constraining Grain Alignment Mechanisms 

The mechanism explaining how non-spherical dust grains align with a magnetic field has 



been sought for ever since Hall (1949) and Hiltner (19491 discovered dust-induced starlight 



extinction polarization. This quest has been highly controversial, and numerous mechanisms 
have been proposed and developed to various degrees over the years (see, e.g., the review by 



Lazarian 2007). Among those is the paramagnetic Davis-Greenstein mechanism that has been 



matured through intensive work ( 


Jones &: Spitzer 1967 


Purcell 1979 Spitzer & McGlynn 


1979 Mathis 


1986 Roberge et al. 




L993 Lazarian, 1997 


Roberge &: Lazarian 1999 ) since its 


introduction by Davis &: Greensteiu 


L (1951). However, this "textbook solution" to the grain 



alignment problem has several shortcomings, among which the fact that thermal flipping 



prevents grains with effective radii a < 1 fim from aligning ( Lazarian 2007 1 , which is hard to 



reconcile with fits to starlight polarization data (Draine Sz Fraisse 2008). Radiative torques 



appear to not only solve this problem, but to also explain why, and not only how, large grains 



(a > 0.1 fxm) align while small grains don't (as first established by Kim &; Martin 1995). 



This promising mechanism was first discovered by Dolginov Sz Mitrofanov (1976). Consid- 
ering a grain exhibiting a difference in optical cross-section for right-handed and left-handed 
photons, they noticed that scattering of unpolarized light by the grain resulted in its spin-up, 
which leads to its alignment with the surrounding magnetic field. However, they did not have 
the necessary tools to quantify this effect and, as a result, their pioneering work was mostly 



ignored in the following 20 years. Draine & Weingartner (19961 realized that this mechanism 



could be studied with a modified version of the DDSCAT code by Draine Sz Flatau (1994). 



In particular, this work, along with follow-up studies by Draine &: Weingartner (1997) and 



Weingartner Sz Draine ( 2003 ) , demonstrated that the magnitude of radiative torques is very 



substantial for the grain shapes they study. It is worth pointing out that the spin-up of grains 



by radiative torques was later observed in laboratory conditions by Abbas et al. (2004). 

However, radiative torques have the inconvenience of depending on many parameters, such 
as grain shape, grain size, radiation wavelength, grain composition, and the angle between the 
radiation direction and the local magnetic field. It is therefore difficult to compute their effect 
in general, and empirical studies have so far been focused on demonstrating the radiative 
torque effect for specific values of these parameters. As a result, it has been claimed that the 
predictive power of this theory was limited, which is true if a "brute force" approach is used 
to explore the relevant parameter space. An alternative is to use analytical models, such as 



the one proposed by Lazarian Sz Hoang ( 2007 ) , to describe the effects of radiative torques 
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6.1 Analytical Model of Radiative Torques 

Lazarian & Hoang (2007) demonstrated that the simple model shown in Fig. |4] can be used to 
accurately calculate the radiative torques exerted by an incoming radiation on a dust grain. 
In this model, grains are assumed to be ellipsoidal, and each grain is equipped with a mirror 
attached to its side. A grain can be both "left-handed" and "right-handed" (from the point of 
view of the incoming radiation) depending on the orientation of its mirror with respect to the 
incident light - the grain shown in Fig. |4]is left-handed, and would become right-handed if the 
mirror were rotated by 90° in a. Radiative torques are then computed under the assumption 
that geometric optics apply, and compared to numerical calculations performed by DDSCAT. 




= 



Figure 4: Lazarian & Hoang (20071 model of a "helical" grain consisting of a splieroidal grain with an 



inclined mirror attached to it. This model can be used to simplify the calculation of radiative 
torques (see text), k indicates the direction of propagation of the incoming radiation, and 
G the angle between the incident light and the polar axis of the ellipsoidal grain. 



The comparison is done for light propagating in the direction of the grain's short axis ei, 
while the grain's angular momentum J is in the (61,62) plane. For the purpose of modeling 
grain alignment, only torques in that plane then matter - a torque in the 63 direction would 
induce grain precession negligible compared to the Larmor precession of the grain in the 



interstellar magnetic field (see Tab. 1 of Lazarian (20071 for the time scales involved). As 
a result, the model only needs to accurately reproduce the components Qei and Qe2 of the 
torque Q along ei and 62, which it does. Defining O as the angle between the magnetic 



field and the radiation direction, Lazarian Sz Hoang (20071 indeed showed that the deviation 



of the torques calculated numerically for irregular grains by DDSCAT from their analytical 
predictions is always modest at the wavelengths of interest. As an illustration. Fig. |5] shows 



(A2)(Qe2) 



7r(Q--)2 



/^model 
^e2 



(©) 



de. 



(5) 



which characterizes the deviation in the 62 direction, as a function of wavelength for various 
grain shapes and sizes. The value of remains small no matter the values of these parameters, 
which indicates a remarkable agreement between analytical model and numerical calculations. 
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In fact, Lazarian Sz Hoang (2007) showed that, for a given grain, radiative torque ahgnment 
is fully determined by the ratio q = Q'^^^/Q^^^ = Qei(0)/(5e2(7r/4). This enormously simpli- 
fies the calculations of radiative torques since one does not need to calculate two functions, 
Qei(0) and Qe2{G))i but only two numbers, Qei(O) and (5e2(7r/4), to characterize the align- 
ment, q is therefore as important to grain alignment as grain axial ratios are to determining 
the polarized radiation by aligned interstellar dust grains. 
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Figure 5: Numerical comparison of the radiative torques calculated with DDSCAT for irregular grains 



with the results provided by the analytical model ("AMO") described in Sec. 6.1 using the 
wavelength-dependent estimator defined by Eq. ([5|. a^^g refers to the radius of a sphere 
with the same volume as the ellipsoid used in the analytical model (the index is omitted in 



the text), and the considered shapes are DDSCAT shapes (see Lazarian & Hoang 20071 



6.2 Alignment Configuration 

Observations indicate that interstellar grains tend to align with long axes perpendicular to the 
local magnetic field, a fact once frequently used to argue that the Davis-Greenstein mechanism 
had to be responsible for the alignment. Radiative torques also explain this configuration. 

Interstellar grains indeed experience internal relaxation that tends to make them rotate 
about their axis of maximal moment of inertia, which, in the case of an ellipsoidal grain, is the 
grain's short axis. Thus, the grain's angular momentum J is perpendicular to its long axis, and 
it is therefore sufficient to follow its dynamics to determine the alignment configuration. Let us 
call the component of the torque parallel (respectively, perpendicular) to J the spin-up torque 
H (respectively, the alignment torque F). The angular momentum J is precessing about the 



magnetic field B due to the magnetic moment of the grain (Dolginov Sz Mitrofanov, 19761. It 



is possible to show that, when the angle ^ between J and B tends to zero (or, by symmetry, vr), 
the alignment torque F averaged over one Larmor precession of J about B also tends to zero 
(see Lazarian &; Hoang 2007 1 . In other words, (F)^ vanish as ^ — > or vr irrespectively of the 
radiative torques' functional form, and the points ^ = and vr are attractors. 

This qualitative argument is quite general, but it does not address the question of whether 
there are other stationary points, and in particular of whether the alignment can also happen 
with J perpendicular to B. To answer this question, one has to consider the behavior of the 



functions Qei{Q) and Qe2(0)- The analysis in Lazarian &; Hoang (20071 shows that there is, 
indeed, a range of angles between the direction of radiation and the magnetic field for which 
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grains tend to align in the "wrong" way, i.e., with long axes parallel to the magnetic field. 
However, this range is very narrow around 7r/2, and therefore not typical of conditions in the 
interstellar medium. Moreover, such an alignment corresponds to positions for which the spin- 
up torque H is negative, which induces grain alignment with low angular momentum. Grain 
thermal wobbling at low- J attractor points leads to variations in ^ that typically exceed the 



range of angles in which "wrong alignment" happens (Lazarian 1994 Lazarian Sz Roberge 



19971 . Therefore, grains always align with long axes perpendicular to the magnetic field. 



6.3 Testing Quantitative Predictions with CMBPol 



In Sec. 6.1 we pointed out that, when induced by radiative torques, grain alignment is fully 
characterized by the value of the parameter q. For typical interstellar conditions, this process 
dominates over all other alignment-inducing torques, and alignment of interstellar dust grains 
should therefore be governed by q independently of the magnitude of radiative torques. 
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Figure 6: Type of available attractor points in the (q, tp) parameter space (see text for the definition 
of these parameters). When both high- J and low- J attractor points are present, grains 
eventually perfectly align in the state of highest J. If only low-J attractor points are 
present, the alignment is partial. "ISRF" refers to the local interstellar radiation field, and 



the various shapes considered here are DDSCAT shapes (see Lazarian & Hoang 20071. 



Fig. |6] shows how a given value of q can be translated into information about grain align- 
ment. In particular, depending on the angle ■0 between the incident radiation and the local 
magnetic field, grains align with ^ taking the value of high- J or low-J attractor points. In 
the former case, alignment with the local magnetic field is perfect. However, when ^ takes the 
value of a low-J attractor point, the alignment is unstable, which leads to a degree of align- 
ment of only ~ 20% to 50%. This situation is actually the most common for relevant values 
of q, meaning that most grains align suhthermally. This conclusion is quite different from the 



assumption by Draine k, Weingartner ( 1996 ) that in the presence of radiative torques most 
interstellar grains should rotate with Trot 3> Jgas, and therefore align suprathermally with B. 
When superparamagnetic grains, i.e., grains with enhanced paramagnetic relaxation, are 



considered, the conclusion is quite different. These grains, invoked by Jones &: Spitzer (1967) 
in the context of paramagnetic alignment (see also arguments in support of their existence 
in 



Bradley 


1994: 


Martin 


1995 



points when subjected to the diffuse interstellar radiation field, as for ordinary paramagnetic 
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grains, but a more stable high- J attractor point also always exists. Therefore, if a grain aligns 
at a low- J attractor point, gaseous collisions ensure that its state will change to that of the 
high- J attractor point. Thus, rather unexpectedly, intensive paramagnetic relaxation changes 
the rotational state of dust grains, enabling them to rotate more rapidly, and leading to perfect 
suprathermal alignment with the local interstellar magnetic field. 

Although we mentioned that radiative torques dominate over all other alignment-inducing 
torques in the diffuse interstellar medium, this might not be the case for paramagnetic grains 



undergoing H2 formation at their surface. The resulting "pinwheel" torques (Purcell, 19791 



can indeed create new high- J attractor points (Hoang Sz Lazarian 20081. As a result, pro- 



vided that H2 torques are as strong as considered in Purcell ( 1979 ) and subsequent papers by 



Spitzer &; McGlynn (19791, Lazarian (1995), and Lazarian Sz Draine (19971, higher degrees 
of polarization are expected in regions of higher atomic hydrogen content. Detection of these 
variations would be highly interesting since it would rule out the superparamagnetic dust 
grains hypothesis for which the alignment should be perfect irrespectively of other factors. 
Thanks to the existence of these quantitative predictions, the theory of grain alignment is 



now becoming testable. Lazarian (2007) reviews some of the first tests the theory of radiative 



torque alignment has already passed. Since then, it has also been tested in dense environments 



such as molecular clouds (Andersson & Potter 2007 Whittet et al. 2008). The current 



agreement between theory and observations is encouraging, and the mere fact that quantitative 
predictions exist that can be tested shows the tremendous progress made in the study of 
grain alignment over the course of the last few years. However, additional tests are required, 
and CMBPol offers a great opportunity to further constrain theories of grain alignment. By 



improving our understanding of the structure of the Galactic magnetic field (Sec. 4.2 ), and our 



ability to measure polarized dust emission (Sec. 5.5) and morphologically trace foregrounds 



(which we discuss in Sec. 8.3), CMBPol will indeed make it possible to test radiative torque 
alignment predictions at multiple frequencies and in various environments. 

In the meantime, observations at smaller scales, and later on with the Planck satellite, 
will improve our knowledge of grain alignment. With a better understanding of the degree of 
alignment of various types of grain, it will become possible to make more robust predictions of 
the expected polarized dust emission at CMB frequencies by combining fits to data at other 



wavelengths with theory, therefore improving on the Draine & Fraisse (2008) results. These 
improved templates will be helpful in generating more accurate sky models used as input to 
study how to best separate foreground emissions from cosmological signal (for more details on 
this topic, we refer the reader to the companion paper by punkley et al.||2008a| . 



7 Heliospheric Dust Emission 

7.1 Basic Properties from In Situ Observations 

Circumheliospheric interstellar material (CHISM) flows through the heliosphere with a helio- 
centric velocity of 26.3 km/s from the upwind direction i ~ 3.5°, b ~ 15.2°. Approximately 
0.3% to 0.9% of the CHISM mass is contained in interstellar dust grains with masses ranging 
from 10-15 to 10-1° 



Griin et al. 


1994 


Frisch et al. 


1999 


Landgraf et al. 


2000 



Slavin &; Frisch 2008), corresponding to radii a between ~ 0.04 //m and 2 //m for spherical 
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grains with density 2.5 g/cm'^. The dynamic interaction of an interstehar dust grain (ISDG) 
with the hehosphere depends on the ratio of radiation pressure to gravitational force (/3), as 
weh as on the Lorentz force, which is a function of grain velocity, charge-to-mass ratio {Q/m), 
and solar wind magnetic field. Consequently, the distribution of ISDGs in the hehosphere 
varies between upwind and downwind, i.e., between Northern and Southern ecliptic regions, 
and with the 22-year solar magnetic activity cycle. Ulysses observed a varying flux of grains 
through the inner hehosphere of 0.6 - 2 x 10~® cm~^ s~^ (Landgraf et al. 2003). The com- 



position and optical properties of ISDGs interacting with the hehosphere are consistent with 
astronomical silicates such as olivine; organic particles may be excluded by the overabundance 
of carbon in CHISM gas ( Landgraf et al. 1999 Slavin Frisch[ 2008 1. However, Helios data 
inside 1 AU showed large ISDGs (radii between 1 and 2 yum), with mixed composition grains 
containing both silicates and carbonaceous material (Altobelli et al. 20061. 




Figure 7: Simulations by Landgraf (20001 of the distribution of interstellar dust grains with radii 
a ~ 0.2 fim in the X-Z plane (see below) are shown for the years 1996, 1998.5, 2001, 
2003.5, 2006, 2008.5 (upper left through lower right). The 1996 figure (upper left) shows 
the grain distribution during solar minimum defocusing conditions, while the 2006 figure 
(lower middle) shows the distribution during solar minimum focusing conditions. ISDGs 
move through the hehosphere at velocities of about 5 AU/yr, so the dust distribution in 
each panel reflects the influence of the solar cycle over several years prior to that date. 
The 1996 solar cycle phase is similar to that expected for a 2018 CMBPol launch date. 
Individual plots are 80 AU x 80 AU, with the Sun at the center. The X-axis is directed 
towards the upwind direction of the CHISM flowing through the hehosphere, while the 
Z-axis is perpendicular to the X-axis and nominally shows the solar rotation axis since the 
upwind direction is 5° from the ecliptic plane. The color scale is given in units of the dust 
density in the diffuse interstellar medium, and ranges from (deep purple) to 4 (red). 



The dependence of grain propagation on Q/m and /? leads to several characteristics of the 



spatial distribution of ISDGs in the hehosphere (Kimura & Mann 1998 Frisch et al. 1999 



Landgraf et al. 2003 Frisch, 2007 Kriiger & Griin 20081: (a) There is a severe reduction 
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in the number of grains with radii a < 0.2 fim penetrating to within 5 AU of the Sun when 
compared to the initial interstellar power law distribution (see, e.g., the simulations shown in 
Fig. (b) Large grains, a > 0.5 //m, on hyperbolic orbits are gravitationally focused into a 



tail downwind of the Sun (Grogan et al. 1996 Landgraf 20001. (c) A gap in the distribution 
of ISDGs with mass between 10~^^ and 10~^^ g is found ~ 2 - 4 AU from the Sun due to the 
effect of radiation pressure (Landgraf et al. 19991. (d) Small grains with masses m < 10^^'^ g 



are alternately focused towards, or away from, the ecliptic plane depending on solar magnetic 
polarity (Grogan et al. , 1996 Landgraf 2000). (e) There is a time lag of ~ 4 - 6 years 
between solar cycle phase and its influence on inner heliosphere ISDGs because of the nominal 
5 AU/yr grain propagation velocity through the heliosphere. (f) The CHISM contains large 
ISDGs, m > 10^^^ g, that are difficult to explain in the context of the extinction properties 



of more heavily obscured sightlines (Frisch et al. 1999 Landgraf et al. 2000 Draine 20081. 
Only properties (b), (e), and (f) above are relatively invariant with the solar cycle phase. 

Several current or proposed spacecraft missions have the potential to return data on ISDGs 
in the heliosphere by the time CMBPol flies. The only spacecraft presently capable of making 
in situ measurements of ISDGs beyond 1 AU from the Sun is Cassini. In situ data collected 
during the same time period by Ulysses and Cassini flnd consistent results, detecting grains 
with radii ~ 0.4 jim on hyperbolic orbits with fluxes ^ 2 x 10~^ cm^^ s~^ (see the analysis 



by Altobelli et al. , 2007). The Cassini mission has now been extended to 2010; the further 
extension of this mission to ~ 2016 (depending on spacecraft consumables) would provide in 
situ data on ISDGs for a similar solar magnetic polarity as expected for the CMBPol launch 
time. Proposed interstellar dust missions include DUst Near Earth (DUNE) ( |Grun Sz Srama 



2006[ ) and SAmple Return of IntersteUar Matter (SARIM) ( [Srama et al.| [2008j). DUNE wiU 
have the capability of acquiring the mass distribution, charge state, elemental and isotopic 
composition, and velocity vectors of ISDGs near 1 AU, while SARIM will return interstellar 
and interplanetary dust grains for laboratory analysis, providing "ground-truth" on ISDGs. 



7.2 Emission Characteristics and CMBPol Observations 

The heating of interstellar dust grains 1-90 AU from the Sun is dominated by solar radiation, 
yielding a grain temperature in K at distance R from the Sun of ^ 278 R^^'^, where R is in 
units of AU. The emission of these grains is expected to be ~ 10"^ of zodiacal light emission 



(Mukai 1981 Frisch et al. 1999). For some inner heliosphere viewing conflgurations and solar 



cycle phases, strong asymmetries in ISDG emissions are expected for leading versus trailing 
sightlines with respect to the Earth's orbital motion, and high versus low ecliptic latitude 
lines-of-sight. Grogan et al. (1996) predict emissions of ~ 0.1 MJy/sr at 12 ^um for sightlines 



that view the hotter ISDGs near the Sun during a focusing phase. Smaller grains at low solar 



elongation angles dominate this emission. Fixsen Sz Dwek (2002) flnd the ratio of the 12 fim 
to 340 /xm flux to be /i2/-^340 ~ 500 for zodiacal light. Full models of the emissions from 
heliospheric interstellar dust are under development (Brahman et al., in preparation), but, 
based on this value of I12/I3AQ, we expect that the emission from interstellar dust grains in 
the heliosphere could be as bright as 0.2 kJy/sr during a focusing cycle. 

The distribution of ISDGs inside of the heliosphere is highly asymmetric around the Sun. 
Figure [7] shows simulations of the distribution of a = 0.2 /im ISDGs in the X-Z plane (see 
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caption for definition) at ~ 2.5 year intervals starting in 1996 (Landgraf 20001. During tlie 



year 2018, which could see the launch of the CMBPol satellite, conditions similar to those 
shown for the year 1996 are expected. For comparison, the COBE data were collected at the 
end of the 1980's focusing cycle, while WMAP, launched in 2001, has witnessed the phases 
shown in the last four panels of Fig [7] (starting by the top right panel). The flux from 0.2 fim 
ISDGs in the inner heliosphere varies by over an order of magnitude between solar minimum 



and solar maximum conditions (Landgraf et al. 2003). As a result, a CMBPol mission is likely 



to detect time variability in the emission from heliospheric dust if it is operated for a long 
enough time. This result would be interesting by itself, since it would probe both the magnetic 
conditions in the heliosphere and the properties of ISDGs, but it would also be useful to the 
CMB analysis by signaling a possible contamination by heliospheric dust emission. 

8 Anomalous Emissions 
8.1 Anomalous Dust Emission 

Prior to the launch of WMAP, detailed intensity maps of the Galaxy were available at 408 MHz 



(Haslam et al. , 19821 and 100 //m (Schlegel et al. 19981. The morphologies of these maps 



were notably different, reflecting the different emission mechanisms that are dominant at those 
frequencies: at 408 MHz, the diffuse Galactic emission is dominated by synchrotron, with large 
characteristic morphological features such as Loop I (the North Polar Spur); at 100 ^m, the 
diffuse emission is generated primarily by dust grains emitting blackbody radiation as they are 
heated by the ambient interstellar radiation field. When these two maps were extrapolated to 
WMAP frequencies, assuming reasonable power law models, it was thought that the 94 GHz 
map would look morphologically very similar to the 100 ^m map since thermal emission from 
interstellar dust dominates the diffuse Galactic signal at these frequencies, while the 23 GHz 
map (away from bright HIT regions and filaments where free-free emission dominates) would 
closely resemble the synchrotron dominated 408 MHz map. However, the WMAP 23 GHz 



map looks much more like the 100 fim map than like the 408 MHz map (Hinshaw et al. 2008 1 . 



The origin of this dust-correlated emission in the low frequency WMAP channels is still 
not completely certain, though there is significant evidence that it arises from the smallest 
interstellar dust grains excited into rotational modes and emitting electric dipole radiation 
dDraine k Lazarian[ p98a||bl ). This " spinning dust" model is currently the leading hypothe- 



sis, although "magnetic dust" emission (Draine & Lazarian 19991 and "dust-correlated syn- 



chrotron" emission (Bennett et al. , 20031 are sometimes offered as alternative explanations. 



Both total intensity and polarization observations can provide useful information as to which 
model can best explain the observed anomalous dust-correlated emission in the WMAP data. 



8.1.1 Sorting through Candidates 

Total Intensity Prior to WMAP^s launch, there was already evidence against the hypoth- 
esis of dust-correlated synchrotron emission. Using the Tenerife 10 and 15 GHz data along 



with the COBE/DMR maps, de Oliveira-Costa et al. (19991 showed that the dust-correlated 



emission does not rise with decreasing frequency as expected for synchrotron, but rather falls 
off below 30 GHz, a feature which is characteristic of spinning and magnetic dust emissions. 
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Furthermore, Finkbeiner et al. (20021 identified at least one confirmed example, L1622, of a 



cloud in which the spectrum of emission from 5 to 10 GHz resembles a spinning dust spectrum. 



As part of the first-year WMAP data release, Bennett et al. (20031 presented maps of the 



spectral index of synchrotron emission derived under the assumption that all of the anoma- 
lous emission is synchrotron. These maps have large diffuse regions in which the spectral 
index of synchrotron is /? > —2.6, whereas, if the synchrotron emitting electrons are accel- 
erated through Fermi acceleration by supernova shocks, the injected spectral index should 
be /3 ~ —2.5 (see, e.g., Gieseler &; Jones 2000), which softens to /3 ~ — 3 when taking into 
account diffusion energy losses (the total energy losses from synchrotron and inverse Compton 
scattering are proportional to E'^). Such a discrepancy seems to invalidate their initial as- 
sumption, a conclusion reinforced by the fact that a spectral index of /3 ~ —3 has been found 



in more rigid fits to the WMAP temperature data (Dobler Sz Finkbeiner 2008a), by combining 



observations at 408 MHz and 1.42 GHz to predict the emission at frequencies above 20 GHz 



(La Porta et al. 2008), and by using only low frequency WMAP polarization data (see Kogut 



et al. 2007 and the Polarization discussion below for more polarization related arguments). 



Targeted observations of individual diffuse clouds have also proven fruitful (see, e.g., iDick- 



inson et al. 2006), and there is evidence that when combining WMAP data with data sets 



at other frequencies, diffuse, inverted spectrum emission (as opposed to the monotonically 



decreasing spectrum expected for synchrotron) is seen (Boughn Sz Pober, 2007). Additionally 



the Green Bank Galactic plane survey (Finkbeiner et al. 2004) presented very strong evidence 



for diffuse inverted spectrum emission between 8.35 and 14.35 GHz outside HII regions. 



Finally, Dobler Sz Finkbeiner (2008b) and Dobler et al. (2008) showed that a spinning 



dust spectrum is completely recoverable within the WMAP frequency range. However, this 
WIM ( "warm ionized medium" ) spinning dust emission correlates more strongly with maps of 
Ha emission (Finkbeiner 2003[ ) than with maps of 100 ^m emission. Since Ha maps represent 



an emission measure, and since spinning dust emission is also a density squared process (Dobler 



et al. 2008), Ha should trace spinning dust emission as well. Magnetic dust emission, on the 



other hand, should be proportional to density and would not be traced by Ha. 



Polarization Perhaps the strongest evidence against the dust-correlated synchrotron hypoth- 



esis is derived from the WMAP polarization maps presented by Page et al. (2007) and Kogut 
et al. (2007). Synchrotron emission is thought to be highly polarized, typically ~ 40 - 70% in 



the presence of reasonably ordered magnetic fields, although the polarization signal can effec- 
tively be eliminated for lines-of-sight going through regions with turbulent magnetic fields (see 



Sec. 4.2 and, e.g., Miville-Deschenes et al. , 2008 ). The polarized fraction of spinning dust emis- 



sion, on the other hand, is expected to be smaller than 3% at the lowest WMAP frequencies 



(Lazarian & Draine 2000). However, the WMAP maps of the polarized emission at 23 GHz 



look morphologically very different from the corresponding total intensity map. In particular, 
the polarization map more closely resembles (but is not perfectly traced by) the 408 MHz map 
and does not correlate very strongly with the 100 fim dust map. The implication is that there 
is a highly polarized component that morphologically looks very much like synchrotron and 
a low polarization component that looks very much like dust. The polarization properties of 



magnetic dust are not well known (see Draine & Lazarian, 1999). 
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8.1.2 Towards a Spectral Identification 

Presently, the study of anomalous dust emission is mostly aimed at targeted searches of promis- 
ing regions (e.g., dust clouds and HII regions) that may exhibit evidence for anomalous emis- 



sion spectrum. These studies are being undertaken in both intensity (Dickinson et al. 20061 



and polarization (Mason et al. 2008). The polarization studies in particular will prove very 



useful in determining just how much of an issue certain emission mechanisms will present 
when attempting to clean polarization maps to get a cosmological signal. However, we stress 
that determining the polarization properties of spinning dust does not exclude the possibility 
that other emission mechanisms (such as magnetic dust) might exist with very significant 



polarization fractions exhibiting a strong wavelength dependence (Draine Sz Lazarian 19991 



With the Ha-correlated spinning dust signal detected in the WMAP data (see Sec. 8.1.1 1 
progress has been made towards measuring specific parameters of the dust grains and envi 



ronment (Dobler et al. 20081, but is limited by the relatively narrow frequency coverage of 



WMAP. Near-term and future missions will address this issue. For example, the C-Band All- 



Sky Survey (Pearson 20071 will map the sky in temperature and polarization at 5 GHz. This 



survey should determine if the dust-correlated emission continues to rise from 23 GHz down 
to 5 GHz (indicating synchrotron) or if it has fallen (indicating spinning or magnetic dust). 

The Planck satellite will now launch shortly and with it will come high resolution maps 
at 545 and 857 GHz. At these frequencies, the foreground emission is totally dominated by 
thermal dust emission and not only will this give insight into the polarization properties of 
dust, but it will also significantly facilitate foreground cleaning, especially if the polarized 
fraction of dust emission is independent of wavelength (see Sec. [s]). Obtaining cleaner maps 
is essential to deriving an unbiased spectrum of foreground emission at anomalous emission 
frequencies, i.e., in the range 20 - 100 GHz (Hinshaw et al. , 2007; Dobler Sz Finkbeiner , 2008a I . 



Additionally, if the Absolute Spectrum Polarimeter (ASP, PI: A. Kogut) is funded, it 
will provide excellent frequency coverage from 30 GHz to 5 THz, with over 300 spectral bins 
dedicated to characterizing foregrounds. Many spectral channels are highly desirable since the 
primary goal in foreground science is to characterize the spectrum of each emission component. 
Furthermore, the more well known the foreground spectra, the easier it will be to clean maps 
of the foreground emission to measure the cosmological signal. The shortcoming of ASP is 
that the angular resolution is a 2° tophat, a resolution considerably lower than even WMAP^s. 



8.2 The WMAP "Haze" 



In an analysis of the WMAP temperature maps, Finkbeiner (2004 1 detected an excess of 



diffuse microwaves coming from an extended region around the Galactic center. This emission 



has also been detected in the three-year temperature data by Dobler & Finkbeiner (2008a I 



and Bottino et al. (20081. When CMB polarization data is not taken into account, this 



excess emission, coined "the haze," exhibits a spectrum too soft to be free-free emission, but 
too hard to be synchrotron from supernova shock accelerated electrons when diffusion losses 
are accounted for (Dobler &; Finkbeiner, 2008a I . It was therefore proposed that the haze 



could be synchrotron emission from a population of hard electrons surrounding the Galactic 



center. However, Gold et al. (20081 do not find any evidence of this hardening in WMAP^s 



observations of the polarized Galactic synchrotron emission (see, e.g., their Fig. 16). Although 
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the WMAP polarization data is not yet sufficient to convincingly rule out the synchrotron 
nature of the haze, this discrepancy is worth investigating. 

As with all inferences of foreground spectra from the WMAP data, there is substantial 
systematic uncertainty in the haze spectrum due to the CMB cross-correlation bias. Any CMB 
estimator is indeed contaminated by some residual foregrounds that then get subtracted off of 
each channel prior to the foreground analysis. As a result, the inferred foreground spectra are 
in turn contaminated by a CMB spectrum. The amplitude of the bias can be hard to estimate 
given that it is related to the chance morphological correlation of the foregrounds with the 
CMB, and, consequently, the haze spectrum is relatively uncertain. Measurements at more 
frequencies will therefore prove crucial to determining the fundamental nature of the haze. 

Thus, the upcoming experiments discussed in the context of spinning dust emission can also 
play a role in characterizing the haze. In particular, with its improved foreground removal abil- 
ities, Planck will lead to a significant decrease in the amplitude of the CMB cross-correlation 
bias. Moreover, the amplitude of the haze in the C-BASS 5 GHz maps will distinguish be- 
tween spinning dust and synchrotron type spectrum for this excess emission. If these two 
surveys confirm that the haze is compatible with hard synchrotron emission, the Fermi satel- 
lite should see 7-rays from inverse Compton scattered starlight (Dobler, in preparation), which 
would confirm the synchrotron nature of the haze. Lastly, the many frequency channels of 
ASP may allow us to trace out the synchrotron spectrum (and thus the energy spectrum of 
the electrons) as a function of position on the sky with a resolution of a couple of degrees. 



8.3 CMBPol as a Morphological Tracer 

The CMBPol satellite has the potential to combine the most promising aspects of upcoming 
experiments and shed new light on both anomalous dust-correlated emission and the haze. 
Although the final design may only have a limited number of frequency bands, the large-scale 
spectral information provided by a survey like the ASP would indeed enable the separation 
of foreground emission and cosmological signal with high accuracy, and in particular, with 
insignificant CMB cross-correlation bias. Thanks to the expected very high spatial resolution 
of CMBPol, one could then trace out the spectrum of individual components morphologi- 
cally, which not only could provide clues regarding the origin of the haze, but is essential to 
constraining the parameters of an environment-dependent emission such as spinning dust. 



9 Foregrounds and CMB Lensing 



Foregrounds are a potential source of systematic error for lensing of the CMB. Quadratic 
lensing reconstruction with CMB polarization indeed generally involves taking products of 



the E and i?-mode polarization signals to compute the integral (Hu & Okamoto 2002 1 

0(1)^ / «;(!', 1-I')i5;(l')i?(l-I')d2l', 



(6) 



where is the estimator for the lensing potential (/i>(l), w a weight function, and 1 the 
wavevector. There are other reconstruction techniques that give smaller error bars (see, e.g.. 
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Hirata &: Seljak, 20031, but, for the purpose of this discussion, we focus on the quadratic 



methods as they are easier to analyze. In general, three types of measurements concern us: 

• the galaxy-lensing potential cross-power spectrum, 

• the lensing potential auto-power spectrum, and 

• the low-^ B-mode power spectrum after correcting for lensing. 

The last of these is the most directly important for the community's quest to detect a pri- 
mordial gravitational wave signal, while the first two are relevant for low-redshift (e.g., dark 



energy and neutrino mass) studies . The companion paper by Smith et al. (2008) focuses on 
the lensing auto-power spectrum. We assess the foreground situation for this particular case. 

As Eq. Q involves a product of E and i?-mode polarizations, the estimator's power 
spectrum is a trispectrum containing two factors of the E-mode polarization and two of the B, 
with both "disconnected" terms (involving the product of two power spectra) and "connected" 
(i.e., non-Gaussian) terms. The disconnected terms depend on the product of the CMB C^^ 
and the foreground power spectrum C^^. The latter are much harder to estimate as they 
depend on the polarized foreground trispectrum, which is mostly unknown at the present time. 

Since lensing reconstruction uses mainly the small-scale polarization, our expectation is 
that extragalactic foregrounds (point sources) will dominate. In this sense we are fortunate, 
in that for point sources with random polarization angles it is possible to write the relevant 
trispectrum in terms of three quantities: (i) the polarized flux-squared-weighted integral of 
the unmasked source counts, I2', (ii) the fourth-power-weighted integral, and (iii) the an- 
gular clustering power spectrum of the point sources, Cf. While the polarized source counts 
at relevant frequencies are presently unknown, we can estimate reasonable ranges for these 
parameters. More importantly, the dependence on a single function raises the hope that de- 
termination of the foreground parameters may be possible by examining the full configuration 
dependence of the trispectrum, a possibility yet to be investigated. This rather optimistic take 
on the foreground trispectrum should however be tempered by some "known unknowns." 

• Point sources have a different spectrum than the CMB, which means that in a broadband 
CMB map they are not truly "pointlike." 

• Extragalactic microwave sources may not have random polarization angles. For example, 
tidal torques may align disk galaxies with each other, resulting in a net alignment of 
magnetic fields. The current data are insufficient to constrain such a contamination. 

• Even in "clean" parts of the sky, there is Galactic emission which may contain non- 
Gaussian small-scale polarized structure with a complicated configuration dependence. 

These problems will have to be continuously re-assessed as we continue with small-scale 
CMB polarization experiments. For example, Planck will measure the small angular scale 
Galactic dust emission and empirically constrain its E/B structure and configuration depen- 
dence, which cannot be accurately predicted from theory or measured with WMAP. We also 
expect that prior to CMBPol, a number of ground-based experiments will make the first char- 
acterization of the lensing S-mode peak, and provide information on the polarization fractions 
and polarized SEDs of point sources at relevant frequencies. Because of the frequency depen- 
dence of the beams, this information will be crucial even if the use of angular information 
(configuration dependence) turns out to be a viable approach to lensing foreground removal. 
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10 Summary and Conclusions 



The main conclusions of this paper are as follows: 



1 . Our current knowledge of the large-scale structure of the Galactic magnetic field comes 



mostly from radio observations of the Galactic synchrotron emission (Sec. 2.1.1) and 



Faraday rotation (Sec. 2.1.2), with the former providing information about the structure 



of the field in the plane of the sky, while the latter encodes information about the line- 



of-sight component of the field. Other techniques (Sec. 2.1.4), such as observations of 



polarized dust emission, starlight polarization, and Zeeman splitting, complete this large- 
scale picture by probing the field on smaller scales, e.g., in molecular clouds, massive 
star forming regions, and low-dust extinction environments. Despite this combination of 



techniques, the large-scale structure of the field is still poorly constrained (Sec. 2.2.1) 



with, in particular, no real agreement as to the number and location of large-scale 
reversals (Fig. [T|. In addition, very little is currently known about both the small-scale 



field (Sec. 2.2.2) and the halo field (Sec. 2.2.3) 



2. Combining multiple probes of the Galactic magnetic field, such as observations of total 
and polarized synchrotron emissions, free-free emission, and Faraday rotation, to simul- 
taneously fit for various parameters describing the magnetized interstellar medium can 
lift degeneracies between them (Sec.jsj), thereby improving constraints on their values. 

3. Knowing the structure and strength of the Galactic magnetic field is a pre-requisite to 



understanding both the total and polarized synchrotron emissions (Sec. 2.1.1) and the 



polarized dust emission (Sec. 5.6). The progress we will make towards characterizing 



the field before CMBPol flies will therefore be of tremendous help in separating Galactic 



foregrounds from cosmological signals. Combining RM synthesis (Sec. 2.1.3) and radio 



polarimetric surveys with wide frequency and spatial coverage will greatly improve our 
knowledge of the field by mapping it in 3 dimensions and enabling the use of RM grids 



(Sec. 4.1 ). In particular, detailed studies of the number and locations of large-scale field 
reversals, the possible inclination of the field with respect to the Galactic plane, the 
vertical structure of the field, and the properties of the turbulent field on a wide range of 
scales will be carried out by the end of the next decade. Data from the Auger experiment 
may additionally allow probing the large-scale structure of the Galactic magnetic field 



without having to rely on a model for the electron density in our galaxy (Sec. 2.1.4) 



4. CMBPol will shed new light on the structure of the Galactic magnetic field through 
higher resolution and sensitivity observations of the total and polarized synchrotron 



emissions from the Galaxy than WMAP and Planck (Sec. 4.2). These observations 
are unaffected by Faraday rotation and therefore complement RM data by probing the 
sky-projected component of field. However, the sensitivities of Planck and WMAP will 
not be high enough for them to accurately determine the characteristics of the polar- 
ized synchrotron emission at intermediate and high Galactic latitudes. With 10 times 
higher sensitivity, CMBPol would be able to constrain the field at these latitudes, and 
in particular its vertical structure and its halo component. Higher resolution will also 
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allow CMBPol to study turbulent structures in the magneto-ionized medium, and to 
potentially trace magnetic fields from the diffuse ISM into small-scale structures. 

5. Models of interstellar dust (Sec. |5.2[) based on observations of the Galactic extinction 



and polarization of starlight (Sec. 5.1) can be used to predict the diffuse polarized far- 



infrared and submillimeter dust mission from our galaxy (Sec. 5.3). However, predictions 
are different for models involving different species of polarization-contributing dust grains 
(Fig.[3|. Although this might make separating polarized dust emission from cosmological 



signal more complicated (Sec. 5.4), it also offers the opportunity of improving our knowl- 



edge of the properties of interstellar dust grains with CMBPol observations (Sec. 5.5). 



6. Starlight polarization data can provide useful Galactic magnetic field templates to trace 



the polarized dust emission (Sec. 5.6), but few data points are available at latitudes 
of interest to CMB studies, typically |6| > 10°. No currently underway or planned 
experiment addresses this issue, since all of them focus either on regions of the Galactic 
plane or on dense small-scale environments. We believe obtaining this data by the time 
CMBPol flies should be a high priority item. 

7. Alignment of dust grains with the Galactic magnetic field is the physical process behind 
both polarized dust emission and starlight polarization. The radiative torque theory of 
grain alignment (Sec. seems to be able to explain why and how dust grains in the 



interstellar medium align the way observations indicate (Sec. 6.2), but the most exciting 
feature of this theory is its ability to make testable predictions (Sec. 6.3) thanks to simple 
models accurately reproducing complex numerical calculations (Sec. 6.1). CMBPol's 
high sensitivity and resolution, combined with starlight polarization information, will 
therefore allow testing the theory of radiative torques by observing the diffuse polarized 



dust emission from the Galaxy in numerous independent patches of the sky (Sec. 5.5). 



8. The distribution of dust grains in the heliosphere, and therefore the corresponding far- 
infrared and submillimeter dust emission, vary cyclically due to the varying magnetic 
conditions originating from the solar activity cycle (Fig. [7]). Although more precise esti- 



mates are needed, CMBPol could detect these fluctuations (Sec. 7.2), thereby constrain- 



ing the magnetic conditions in the heliosphere and the properties of heliospheric dust. 
9. Two anomalous emissions were detected in the WMAP data. Dust- correlated emission in 



its low frequency channels is likely to be spinning or magnetic dust emission (Sec. 8.1 ) 



whereas the origin of an excess emission of microwaves towards the Galactic center. 



dubbed "the haze" (Sec. 8.2), is still unknown but could be synchrotron radiation from 
a population of hard electrons. A definitive identification of these two emissions likely 
requires observing their spectra, an effort ground-based experiments, e.g., C-BASS, and 



dedicated satellite missions, such as the ASP, will contribute to (Sec. 8.1.2). CMBPol's 
high resolution will enable us to trace these spectra morphologically (Sec. 8.3). This 
is essential to determining the parameters of spinning dust emission, which strongly 
depend on environment, and it might give us clues as to the possible origin of the haze. 

10. Finally, extra-Galactic sources are a potential contaminant for studies of CMB lensing 
(Sec. [9]), but not much is currently known about the polarized source counts at CMB 
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frequencies, or the spectra and intrinsic polarizations of these sources. The input of 
ground-based experiments on these questions wih be essential to improve our measure- 



ments of CMB lensing (see more details in the companion paper by Smith et al. 2008 1 



Detecting a primordial B-mode polarization signal would have dramatic consequences for 
our understanding of the phenomena that took place in the early Universe and led to the 
necessary conditions for it to evolve into the Universe we observe today (see the companion 



paper by Baumann et al. 2008a). Being able to do so will not only require advanced CMB data 



analysis techniques (see the companion paper by Dunkley et al. , 2008a), but also an in-depth 



understanding of how Galactic [and, for some purposes (see, e.g., the companion paper by 



Smith et al. 20081, extra-Galactic] foregrounds behave. Taking this synergy to the next level 



necessitates coordinated and collaborative efforts by both the CMB and the ISM communities, 
each of which, as we hope to have shown throughout this document, will highly benefit from. 
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